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We investigate the properties of membranes under tension by Monte-Carlo simulations of a generic 
coarse-grained model for lipid bilayers. We give a comprising overview of the behavior of several 
membrane characteristics, such as the area per lipid, the monolayer overlap, the nematic order, and 
pressure profiles. Both the low-temperature regime, where the membranes are in a gel Lg/ phase, 
and the high-temperature regime, where they are in the fluid L a phase, are considered. In the 
Lpi state, the membrane is hardly influenced by tension. In the fluid state, high tensions lead to 
structural changes in the membrane, which result in different compressibility regimes. The ripple 
state Ppt, which is found at tension zero in the transition regime between L a and Lg/, disappears 
under tension and gives way to an interdigitated phase. We also study the membrane fluctuations 
in the fluid phase. In the low tension regime the data can be fitted nicely to a suitably extended 
elastic theory. At higher tensions the elastic fit consistently underestimates the strength of long- 
wavelength fluctuations. Finally, we investigate the influence of tension on the effective interaction 
between simple transmembrane inclusions and show that tension can be used to tune the hydrophobic 
mismatch interaction between membrane proteins. 



I. INTRODUCTION 

Biological membranes are made of lipid bilayers with 
incorporated proteins. These barriers define the inside 
and the outside of a cell, separate the functional com- 
partments in cells, and are indispensable for life The 
microscopic surface tension of membranes is usually small 
or vanishes altogether but there may also be situ- 
ations, where membranes are under considerable stress 
due to osmotic pressure differences. For example, ep- 
ithelial cells exposed to transmembrane osmotic gradients 
can be expected to develop a state of tension under phys- 
iological conditions Q. Similarly, osmotically induced 
tension may play a decisive role during conformational 
changes, fission or fusion of cells 4, 5|. Another situation 
where membranes experience stress is under the influence 
of ultrasonic pulses [6| . Applied perpendicular to a lipid 
membrane, shock pulses can promote significant struc- 
tural changes similar to those induced by lateral tension. 
The effect of such pulses on membranes is of considerable 
medical interest. In this context Koshiyama et al. have 
studied phospholipid bilayers under the action of a shock 
wave in atomistic molecular dynamics simulations [?|]- 

Despite the advances in computer technology through- 
out the past decades, atomistic modeling of lipid bilayers 
on length scales of a few nanometer still requires huge 
computing resources or even goes beyond the current ca- 
pabilities of high performance architectures. This mo- 
tivates the use of coarse-grained models. They can give 
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fundamental insights into the physics of a certain system, 
and, moreover, they provide powerful tools for the inter- 
pretation of the behavior of complex systems, like lipid 
membranes 

The aim of the work presented here is to study the 
effect of an externally applied tension on the physical 
properties of a model bilayer, and on the behavior of in- 
corporated model proteins. We employ a generic coarse- 
grained model of amphiphiles developed in a top-down 
approach[l3j. For tensionless systems this model has al- 
ready been used very successfully to reproduce various 
bilayer phases including the symmetric and asymmetric 
" rippled" Pg< states fill . [l5| and to study membrane- 
protein interactions [la ]. Recent simulations on mem- 
branes under mechanical stress have often dealt with the 
formation, structure and stability of hydrophilic pores 
[T3-0 or w hh the influence of surface tension on trans- 
membrane channel stability and function [2(J [H| ■ In this 
paper, we will primarily be concerned with the structural 
changes of pure membranes in response to lateral stresses, 
focussing on unporatcd systems. Since our model ex- 
hibits a rather realistic phase behavior of the model mem- 
brane at different temperatures, we can study different 
membrane states, i.e., the liquid, the ripple, and the gel 
state 

Our paper is organized as follows: First we describe the 
underlying lipid model and outline the simulation meth- 
ods. Then the simulation results are presented starting 
with a phenomenological introduction, where the effects 
of an external tension on the model bilayer in different 
phases are described. Thereafter a quantitative analysis 
of these bilayers is performed and the characteristics of 
the bilayers are examined with respect to their behavior 
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under external tension. In the last part we show the in- 
fluence of tension on the effective interaction potential of 
two simple model proteins. A brief summary concludes 
our paper. 



II. MODEL AND METHODS 

In our model [22j the lipid molecules are represented by 
chains made of one head bead and six tail beads, and we 
have additional solvent beads. Hence we have three types 
of beads, h,t, and s for head, tail, and solvent beads. 
Within the lipid chains, neighbor beads at a distance r 
interact via a finite extensible non-linear elastic (FENE) 
potential 



VfeneM = --e FENE (Ar, 1 
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with the spring constant efene, the equilibrium bond 
length ro, and the cutoff Ar max . The angles 9 between 
subsequent bonds in the lipid are subject to a stiffness 
potential 



y B A(0)=eBA(l-cos0), 



(2) 



with the stiffness parameter £ba- Beads of type i and j 
which arc not direct next neighbors in a chain interact 
via a truncated and shifted Lennard- Jones potential, 







Vc,ij) if r <C Tc.ij 

otherwise. 
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The offset V Ct ij is chosen such that V\,j{r / a.ij) is contin- 
uous at the cutoff r c ^-. The parameter er.y = (<7j +aj)/2 
is the arithmetic mean of the diameters Oi of the inter- 
action partners, and r c jj = 1 Uij for all partners (ij) 
except (tt) and (ss): r cM = 2er tt and r c tSS = 0. Hence 
tail beads attract one another, all other interactions are 
repulsive, and solvent beads do not interact at all with 
each other. This way of modeling the solvent environ- 
ment (the so-called 'phantom solvent') [f3| combines the 
advantages of explicit and implicit solvent models: Like 
implicit solvent models, the solvent environment does not 
develop any artificial internal structure, and it is very 
cheap (in Monte-Carlo or Brownian dynamics simula- 
tions, less than 10 % of the total computing time is spent 
on the solvent beads). Like explicit solvent models, the 
model can be used to study solvent-mediated phenomena 
such as the effect of hydrodynamic interactions on mem- 
brane dynamics. This is much more difficult with solvent- 
free models. We should note, however, that we mainly 
consider static membrane properties in the present work, 
using Monte-Carlo simulations. 

We use the model parameters [22j (Jh = 1 • 1 ft , ?"o = 
0.7 a t , Ar max = 0.2 a t , £ F ene = 1 00 e/of, and e B A = 
4.7 e. At the pressure P = 2.0 e/of, the model repro- 
duces the main phases of phospholipids, i.e., a high- 
temperature fluid L a phase at temperature fc^T > 



ksTm ~ 1.2 e and a low-temperature tilted gel (Lb') 
with an intermediate modulated ripple (Pp>) phase [14j |. 
The energy and length scales can be mapped to SI- 
units [HI by matching the bilayer thickness or, alter- 
natively, the area per lipid and the temperature of the 
main transition to those of DPPC, giving I a t ~ 6 A and 
1 e ~ 0.36 x 1 0~ 20 J. The elastic properties of the mem- 
branes in the fluid state were then also found to be com- 
parable to those of DPPC membranes (l6j . 

Unless stated otherwise, the systems were studied by 
Monte-Carlo simulations at constant normal pressure 
Pn = 2.0 e/of and constant temperature T with peri- 
odic boundary conditions in a simulation box of variable 
shape and size. We put the tension into effect by an ad- 
ditional energy term — TA to the Hamiltonian of the sys- 
tem, where A is the projected area of the bilayer onto the 
xy plane. This alters the lateral components of the pres- 
sure tensor within the membrane. The non-interacting 
solvent particles, which probe the free volume and force 
the lipids to self-assemble, are not affected by this addi- 
tional energy contribution. They ensure that the normal 
pressure P/v is kept fixed at the required value. Thus, we 
are performing Monte-Carlo simulations in an NP^T T 
ensemble with effective Hamiltonian 



H cS = E + P N V - TA - Nk B Tln(V/V ) 



(4) 



where E is the interaction energy, V the volume of the 
simulation box, Vq an arbitrary reference volume and N 
the total number of beads. In contrast to the experimen- 
tal situation, where the lateral pressure of a lipid bilayer 
cannot be controlled very easily (23[, our implementa- 
tion of the tension is straightforward. Since we are in 
full control of the lateral pressure, we can gain insight 
into states and structures of lipid bilayers by means of 
computer simulations, which are difficult to investigate 
in experimental setups. 

In practice, two main types of Monte-Carlo moves were 
proposed and then accepted or rejected according to a 
Metropolis criterion, namely 1) translational local moves 
of the beads, and 2) global moves which change the vol- 
ume of the simulation box or its shape [22T ]. During one 
Monte-Carlo step (MCS) there is on average one attempt 
to move each bead. Since the global moves require rescal- 
ing of all particle coordinates, which is rather expensive 
from a computational point of view, they are performed 
only every 50th MCS on average. The volume and shear 
moves are necessary to maintain the desired surface ten- 
sion. An advantage of our scheme is that the correct area 
per lipid required by the external tension is adopted by 
construction. There is no need for searching the required 
state by testing various values of fixed area per lipid. We 
also implemented flip-flop moves with an inversion of the 
end-to-end vector of a complete lipid chain in its last tail 
bead, but found the acceptance rate of this 'molecular' 
move due to the density of the bilayer far too small. 

In some cases it turned out to be more convenient to 
keep the box height L z fixed and vary only the planar ex- 
tension given by L x and L y . In this case, the number of 
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phantom solvent particles was allowed to fluctuate, i.e., 
additional Monte-Carlo moves were attempted where sol- 
vent particles were removed from the system or randomly 
inserted (semi- grand canonical simulations). The solvent 
chemical potential was set to 

f i s = k B \n(V Q P N /k B T). (5) 

Now, the semi-grand canonical Hamiltonian reads 

H s ,gc = E + P N V-TA-fi s N s -k B T\n [(V/V ) N /N s \)] , 

(6) 

where N s is the momentary number of solvent particles 
in the system. 

When the system sizes were very large (Sec. IIII Dl) . 
the simulations were carried out on parallel proces- 
sors using a domain decomposition scheme described in 
Ref. HH, which ensures that every Monte- Carlo move ful- 
fills detailed balance exactly (in most other decomposi- 
tion schemes proposed in the literature, detailed balance 
is violated at the boundaries between the domains). We 
have checked by comparison with single-processor simu- 
lations that the results were not affected by the paral- 
lclization. 

In this paper the results presented in Section IIII El 
were obtained using the semi-grand canonical solvent 
model, and the diffusivity measurements (Sec. IIII CP in 
the iWT-ensemblc (in the latter case, only local bead 
translations were attempted and all simulations were run 
on a single processor). The other simulations used the 
NP N T T-ensemble or. 

The simulations of stressed membranes were carried 
out for a duration of at least 4 x 10 6 MCS, and we checked 
that the observed quantities did not show any tendency 
for drift. In particular, we have examined in detail the 
long- wavelength undulations of the membranes, which 
are presumably the slowest relaxation modes in the sys- 
tem. After an equilibration time of 2 x 10 6 MCS the re- 
sults from two successive runs of length 2 x 10 6 MCS were 
identical within the error. Hence the analyzed systems 
can be assumed to be in a stable state. 



III. SIMULATION RESULTS 

We have applied tensions of up to T = 3.0 e/of, cor- 
responding to 30mN/m, which is close to the threshold 
for rupture in the fluid phase in our model: The mem- 
branes remained stable up to temperatures k B T — 1.3 e, 
but ruptured at k B T = 1.4 e (see Fig. [T]). This is con- 
sistent with atomistic simulations of Leontiadou et al, 
who found a critical tension for stable porated bilaycrs 
of m 38mN/m in the context of pore formation with an 
atomistic model of DPPC [l8||. In contrast, the bilayer in 
the gel phase can sustain much higher tensions and stays 
stable for values of V up to 4.0e/of and higher. The 
areal expansion and reduction in thickness of the bilayer 
is substantial. In our simulations, we observe an increase 
in area per lipid of more than 40 % and a decrease in 
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FIG. 1: Snapshots (top view) of a rupture event at ksT = 
1.4 e and F = 3.0e/cr t 2 . The underlying curve shows the cor- 
responding evolution of the area per lipid. At about 4.3 erf, 
the area per lipid levels off, because a restriction was imposed 
on the minimal size of the simulation box in the z direction. 
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FIG. 2: Snapshots of bilayer configurations in the fluid phase 
at feT = 1.3 e. Gray scale coding: Light gray molecules 
point upwards from head to tail, dark gray molecules point 
downwards. At increasing tension the two monolayers partly 
interdigitate. The size of the beads are not to scale. 



bilayer thickness of more than 30 % for high tensions, 
where the systems can still be found in non-rupturing 
configurations. In other simulation studies a comparable 
or even larger increase in the area per lipid was observed 
under tension without rupturing of the bilayer. Groot 
and Rabone report an areal expansibility of more than 
70% for mixed membranes before these are finally ripped 
apart [24J. Grafmuller et al. see an areal gain of 60% 
under tension for one of their parameter sets [25j | . 

We begin with giving a qualitative overview over the 
behavior observed in the different bilayer phases. Figs.[5J 
[3l and [4] show configuration snapshots of bilayers under 
tension at the temperatures k B T = 1.3 e, corresponding 
to a fluid state, k B T = 1.2 e, corresponding to a rip- 
pled state, and k B T = 1.1 e, corresponding to a tilted gel 
state. 

In the fluid phase (Fig. [2]), the membrane stretches 
under tension, and the two monolayer leaflets become 
less well separated (right). Thus high tensions change 
the structure of the membrane. This effect is even more 
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FIG. 3: Snapshots of bilayer configurations in the ripple phase 
at ksT — 1.2 e (gray scale coding as in Fig. [2]). On the left 
the tensionless system with a pronounced ripple is depicted. 
This ripple vanishes as tension is applied, giving way to a 
completely interdigitated structure (middle and right figure). 
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FIG. 4: Snapshots of bilayer configurations in the gel phase 
ksT = 1.1 e (gray scale coding as in Fig. [2]). No interdigita- 
tion at increasing tension is observed. 



pronounced in the rippled state (Fig. [3J here the orig- 
inal rippled structure was obtained by cooling tension- 
less equilibrated configurations from the fluid phase at 
fcfiT = 1.3 e down to fc^T = 1.2 e). Under tension, the 
ripple unravels and gives way to an interdigitated phase. 
In contrast, the gel state (Fig. [2} is hardly affected un- 
der tension. The two monolayers remain well-separated. 
Only the average lipid tilt away from the bilayer normal 
is slightly enhanced from (9) = 23° in the tensionless 
state to (0) = 30° at the tension V = 3.0 e/of (snapshot 
not shown). 

After this qualitative overview, we turn to a more 
quantitative analysis of the behavior of bilayers under 
tension. 



A. Global characteristics of pure membranes 

We first consider the area per lipid, which is obtained 
by dividing the projected area of the bilayer in the xy 
plane by half the number of lipids in the system. We 
note that the difference between the number of lipids 
in the upper and the lower monolayer was always very 
small. "Flip- flop" moves are practically never observed 
for systems in the gel phase. In the fluid phase, about 
10 % of the lipids were exchanged between the monolayers 
during 10 6 MCS, but fluctuations of the average number 
of lipids in each monolayer are still less than 2 %. 

We do not observe any dependence of the area per lipid 
on the system size. For the temperature fcgT = 1.3 e, 



TABLE I: Mean area compressibility modulus k A over the 
whole range of tensions. 



we have compared data from four different system sizes, 
ranging from 200 lipids to 7, 200 lipids; the results were 
identical within the statistical error (w 1%). These re- 
sults are in agreement with the findings of Marrink and 
Mark in atomistic simulations um , or with those of Kra- 



nenburg et al. [27J , who studied a coarse-grained model of 
amphiphilic surfactants by a combined DPD and Montc- 
Carlo scheme, imposing the surface tension in a way sim- 
ilar to ours. 

Assuming that the area per lipid A depends linearly 
on the applied tension T, one can calculate the mean 
area compressibility modulus kA using the relation fc^ = 
AoAr/A^4, where Aq is the equilibrium area of the ten- 
sionless membrane. This yields the values of kA listed 
in Table [T] Here only data from configurations which re- 
mained stable for long simulation runs (up to 10 7 MCS) 
have been taken into consideration, i.e., the data for the 
state point ksT = 1.4 e, T — 3.0e/of , which lies beyond 
the rupture threshold, were omitted. 

For membranes in the gel phase (fcsT = 1.0 e and 
ksT = 1.1 e), the number kA fully characterizes the be- 
havior of the area per lipid over the whole investigated 
range of tensions (data not shown). The extensibility of 
the bilayer in the gel phase is significantly smaller than 
in the fluid phase, and constant over the whole range 
of tensions under investigation. The most noticeable ef- 
fect of the tension is the increase in tilt angle of the lipids 
mentioned earlier, which leads to a slightly reduced mem- 
brane thickness. 

The behavior of membranes in the fluid or ripple state 
is more complicated. Fig.[5jshows the corresponding data 
for the area per lipid as a function of the applied tension. 

We first discuss the situation at the temperature 
ksT = 1.2 e, where the tensionless membrane is in 
the ripple phase. As we have already discussed earlier 
(Fig. [3]), the membrane undergoes a phase transition to 
an interdigitated phase under tension. Fig. [51 top, shows 
that this is associated with a pronounced change of the 
compressibility. Up to values of about 1.0 e/ of, the area 
per lipid increases steeply, the resulting value of kA is 
even lower than that of the fluid phase. At high ten- 
sions, r = 1.9e/of and above, the approximate areal 
compressibility is strongly enhanced and comparable to 
values obtained for systems in the gel phase. Thus the 
extensibility of the bilayer changes from fluid-like to gel- 
like under tension. 

In the fluid phase (Fig.[5j bottom), the tension-induced 
structural changes in the membranes are less dramatic, 
but they can still be associated with compressibility 
changes. The data shown in Fig. [5J bottom, suggest a 
subdivision into three different compressibility regimes: 
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FIG. 6: Average overlap parameter (see text for definition) 
vs. tension. Negative values indicate spatially well-separated 
monolayers. For interdigitated bilayers the overlap parameter 
is positive. 



FIG. 5: Area per lipid A vs. tension F for two different tem- 
peratures as indicated, together with linear fits to different 
compressibility regimes. 

Top: As the tension on a patch of bilayer in the ripple-phase 
increases, the areal compressibility changes from fluid-like to 
gel-like. 

Bottom: In the fluid phase, three different compressibility 
regimes are observed, F = 0-1 e/crf , F = 1-2 e/of , and F = 2- 
3e/af . The inset focuses on the tensionless membrane (see 
text for explanation) 



Under tension, the membrane switches from a less com- 
pressible low-tension state to a more compressible high- 
tension state via a highly compressible intermediate. 
Thus the structural changes in the membrane, which were 
observed in Fig. seem to be related to a crossover be- 
tween different membrane states and possibly even reflect 
the vicinity of a hidden phase transition. 

The inset in Fig. [51 bottom, focuses on the limit of 
very small tension/compression. In recent work, kA has 
been determined for the case T = 0.0 e/of with an alter- 
native method, i.e., the detailed analysis of the fluctua- 
tion spectrum of tensionless membranes [l6| . Here, we 
have extracted the areal compressibility modulus of the 
fluid membrane close to tension zero by both compress- 
ing and extending the system slightly within the range of 
r = —0.2 to 0.2e/of . The resulting value for fc^ divided 
by the square of the mean tensionless monolayer thick- 
ness to ~ 3 a t agrees with the value A^/ijj = l.l±0.2e/of 
obtained independently from the fluctuation analysis (sec 
Table |n|. 

To characterize the global structure of the membrane 
in more detail, we next consider the intcrdigitation of 
the two monolayers. It can be characterized in terms of 



the 'overlap parameter' O c hain = (2(Z Z — to)/l z ), origi- 
nally introduced by Kranenburg et al. [27j . The results 
are shown in Fig. [6] They quantify the behavior which 
was already apparent from the snapshots shown earlier 
(Figs. [2]|4]). In the gel phase, where we did not notice 
any overlap of the monolayers even under tension O c hain 
is always negative. In the fluid phase, it becomes positive 
at high tensions. 

Interestingly, this leads to a non-monotonic behavior 
of the orientational chain order in the fluid phase, i.e., 
the nematic order parameter S z = l/2((l z /l) 2 — 1). Here, 
l z is again the z component of the end-to-end vector of a 
lipid chain, and I its full length. As shown in Fig. [71 the 
nematic order first decreases in the low-tension regime 
T < l.Oe/of. At T - 1.5 - 2.0e/cr t 2 , a rather steep 
increase ensues, followed by a plateau in the high-tension 
regime T > 2.0e/af. Under tension, the lipids thus first 
disorder and tilt away from the bilayer normal, which 
leads to an unfavorable packing in the hydrophobic bulk 
of the membrane. As interdigitation sets in, the lipids 
relax and assume once again their preferred order. 



B. Pressure Profiles 

After having discussed these global properties of the 
membranes, we now investigate the effect of the applied 
external stress on the internal stress distribution inside 
the membrane. Stress distribution profiles influence, e.g., 
the permeability of membranes with respect to small 
molecules. To study them, we have recorded the interfa- 
cial tension (or negative stress) profiles 

1 (z) = P zz {z)- l -(P xx {z) + P yy {z)) (7) 
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FIG. 7: Nematic order parameter of chains S z vs. tension. 



in small systems of 200 lipids. The pressure tensor P a p 
is obtained using the virial theorem, 



a/3 



Nk B T 




/3 



(8) 



Here, is the position of particle i, Fj the force acting on 
this particle, AT the number of particles, T the temper- 
ature, and V the volume. The local distribution of the 
pressure along the bilayer normal was obtained by divid- 
ing the system into 50 vertical slabs and distributing the 
pressure contributions onto these slabs according to the 
convention of Irving and Kirkwood [28| . The pressure 
profiles can also be used to cross-check the consistency 
of our approach in Eq. 21 since the integral 



To 



dz 7(2) 



(9) 



has to match the externally applied tension. We have 
checked that this was the case in all simulations. 

Pressure or stress profiles have been reported in various 
other studies @, [M, HI for both mcsoscopic and atom- 
istic simulations. We will briefly summarize the charac- 
teristic features of the total stress profiles and analyze 
their change under tension in our model. The behav- 
ior of the pressure profiles can be attributed to different 
interactions: 

The first positive peak (insets in Figures [5] and [5]) 
arises due to the purely repulsive interactions of the head 
and the solvent beads. Here a zone depleted from (sol- 
vent) beads is formed and the head beads are effectively 
squeezed together in the lateral direction. The first nega- 
tive peak in the head region indicates that the head-head 
interactions are purely repulsive, pushing the system lat- 
erally apart in this section of the bilayer. The second 
positive peak, located in the tail region, results from the 
interplay of the attractive Lennard- Jones interaction be- 
tween the tail beads and the contribution of bond-angle 



and bond-length potentials. Interestingly, the contribu- 
tion of the (attractive) Lennard- Jones (LJ) interactions 
between tail beads is negative: The net attraction be- 
tween tail beads is stronger in the normal direction than 
in the lateral direction. This is compensated and out- 
balanced by the contributions of the bending of the lipid 
segments and the stretching of bonds, which lead to pos- 
itive tension in total. Thus the stress in the hydropho- 
bic portion of the membrane is mostly sustained by in- 
trachain interactions and not, as one might expect, by 
the attractive Lennard-Jones interactions. The negative 
peak in the midplane of the bilayer originates from the 
absence of intrachain interactions in this region, thus the 
effect of the Lennard- Jones interactions takes over and 
the monolayers are effectively glued to each other. 

Under the influence of an external tension, we observe 
a narrowing and shift of the whole profiles, in qualita- 
tive agreement with previous atomistic simulations by 
Gullingsrud and Schulten [2(|. In the gel phase, which is 
already exposed to very high internal stress, the relative 
effect of the external tension is small. The inspection 
of the different contributions to the local pressure shows 
that the external tension leads to a decrease both of the 
Lennard-Jones (LJ) contribution and the bond length 
(BL) contribution. The reduction of the LJ-contribution 
is higher, leading to the observed shift. Due to the stiff- 
ness of the lipids in the gel phase, the bond-angle contri- 
bution to the pressure profiles remains practically unal- 
tered. The change in the overall structure of the pressure 
profiles is also only small. The slight change in membrane 
thickness results in a shift of the outer peaks of the pro- 
files towards the midplane (z = 0). 

The effect of external tension is considerably more dra- 
matic in the fluid phase at fc^T = 1.3 e. Although the 
absolute peak values of the pressure profiles are reduced 
by a factor of about 8 to 10, the relative shifts are much 
more pronounced. First, the shift of the outermost peaks 
towards the midplane is much higher, due to the fact that 
the membrane thickness decreases more strongly under 
tension in this phase. Second, the shapes of the individ- 
ual contributions to the total profile change qualitatively 
under stress, reflecting the structural change from a well- 
separated bilayer to an intcrdigitated structure. At high 
r, the individual terms have a simple structure with sin- 
gle, positive or negative, peaks at the center of the mem- 
brane. Nevertheless, the total stress profile still exhibits 
the oscillatory features described above. 



C. Correlation functions and diffusion of lipids 

Next we consider the effect of the tension on the lateral 
structure of the bilayers and on the dynamic mobility of 
individual lipids. Fig. [TU] shows the two-dimensional ra- 
dial distribution function g m {r) of the (x, y)-coordinates 
of the head groups in the gel and the fluid state for differ- 
ent applied tensions T. In the gel state, g m (f) exhibits a 
series of pronounced peaks, reflecting a high degree of or- 
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FIG. 8: Contributions of the Lennard- Jones (LJ), bond-angle 
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profiles in the gel phase, for the tensionless case (thick lines) 
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the total stress profile, i.e., the sum of all contributions, for 
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FIG. 9: Contributions of the Lennard- Jones (LJ), bond-angle 
(BA), and bond- length (BL) interactions to the total stress 
profiles in the fluid phase, in the tensionless case (left), and 
for the tension T = 2.0 e/of (right). The inset shows the 
total stress profile for the tensionless bilayer (thick line) and 
for T = 2.0e/a t 2 (thin line). 



der (Fig.[TUl left). Under tension, the higher order peaks 
shift to slightly larger distances, reflecting the enhance- 
ment of the area per lipid. The fluid membrane is much 
less structured. The head-head radial correlation func- 
tions within one monolayer decay rapidly already at zero 
tension, and all higher order peaks disappear at higher 
tension (Fig. [TU1 right). In sum, the influence of the ten- 
sion on the lateral structure of the membranes is found 
to be largely negligible. 

The situation is different for the diffusion coefficient 
of lipids. Although Montc-Carlo simulations do not pro- 
vide an intrinsic timescale, one can still obtain valuable 
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FIG. 10: Two-dimensional radial distribution function g m (r) 
of the head groups within one monolayer at ksT — 1.1 e (left) 
and at fcsT = 1.3 e (right). 



information on the diffusional behavior of the lipids from 
NVT Monte-Carlo simulations that employ only local 
bead moves. In our diffusivity studies, the initial con- 
figurations were taken from systems equilibrated in the 
NPnT T-ensemble, but we measured the diffusivity of 
lipids in the iWT-ensemble, i.e., no volume or shear 
moves were carried out during the simulation. We have 
monitored the pressure tensor for the duration of the dif- 
fusion measurements to check that it stayed constant dur- 
ing the simulation. 

In the following, the basic "time unit" is one Monte- 
Carlo step (MCS, the Monte-Carlo 'time scale'), and the 
"time" t counts the number of MCS since the start of 
the simulation. We consider the lateral diffusion constant 
within the membrane (or, more specifically, the projec- 
tion of the membrane into the (xy plane)), defined as [29| 



Dx V (r) = 



6Nt 



(10) 



Here the sum i runs over all lipid heads, rf denotes the 
position of a lipid head i in the xy plane with respect to 
the center of mass of the bilayer, and Ar[ xv ^ its difference 
from one MCS time step to the next without the offset 
that has to be added if the head crosses the periodic 
boundaries. 

In the gel phase (data not shown), no diffusion was 
observed over the whole length of the simulation (r max = 
1.2 x 10 7 MCS). The lipids basically fluctuate around 
their average positions. The width of the fluctuations, 

J ( r (^) 2 ) - (r^y)) 2 ~ 0.04 a t , docs not depend on the 
applied tension. The data for the fluid phase are shown 
in Fig. [TTJ Here the lipids diffuse freely, and the lat- 
eral diffusion constant D xy increases significantly if one 
applies moderate tensions up to T = 2.0e/of (Fig. ITT|) . 
Interestingly, it does not increase further in the high- 
tension regime, beyond T = 2.0 e/of, even though the 
area per lipid and the chain overlap parameter O c ham are 
not yet saturated. 
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FIG. 11: Lateral diffusion coefficient D xy for different tensions 
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It should be noted that the in-plane diffusion constant 
discussed here differs from the true diffusion constant in 
the membrane, due to the presence of membrane undula- 
tions (see also Sec. IIIIDj) . The thermally induced buck- 
ling of the membrane can lead to a substantial out-of- 
plane component to the lipids' diffusional motion, which 
is not captured in our definition of D xy in Eg. 1101 Various 
theoretical studies have addressed this problem [3ll - l33l ]. 
They commonly conclude that membrane fluctuations 
lower the measured in-plane diffusion coefficient. As the 
tension increases, the fluctuations of the membrane are 
suppressed (see Sec. IIIIDj) . Therefore, one may speculate 
that the increase of the apparent diffusivity under tension 
measured in our simulations is, at least in part, caused 
by the reduction of thermal bilayer undulations. How- 
ever, this is not sufficient to fully explain the 20%-cffcct 
observed in Fig. [TTI 

D. Fluctuation spectra 

To conclude the analysis of pure membranes, we study 
the thermal membrane fluctuations, which were already 
mentioned in the previous section. Theoretical consider- 
ations [HJ and experiments [35| have shown that lateral 
tension on fluid bilayers leads to suppression of thermal 
fluctuations, which in turn decreases steric repulsion of 
vesicles and changes adhesive properties. Moreover, the 
Fourier spectra of the height and thickness fluctuations 
provide information on the elastic properties of the mem- 
branes. Therefore, it is instructive to look at the develop- 
ment of undulation, peristaltic and protrusion properties 
of our model membrane under tension. 

To analyze our data, we use an extension of an elastic 
theory by Brannigan and Brown [36j , which we have al- 
ready applied with success to the case of the tensionlcss 
membrane (l6| . We shall not rederive the theory here, 
but merely sketch the main assumptions. Brannigan and 



Brown describe a planar fluctuating membrane as a sys- 
tem of two coupled monolayers surfaces hi^(x, y), of with 
each is characterized by two independent fields, -21,2(2;, y) 
and Ai j 2(a;, y) (h = 2 + A), one accounting for slow 'bend- 
ing modes' and one for fast 'protrusion modes' [13 ■ They 
furthermore make a number of approximations, which 
amount to the assumptions that 

(i) the protrusion modes of the two monolayers are 
independent degrees of freedom, 

(ii) the bending modes 2:1,2(2;,?/) can be rewritten in 
terms of their sum and their difference (z± = 
(2i±22)/2), corresponding to (bending) height and 
(bending) thickness modes of the membrane, which 
in turn decouple. The (bending) thickness modes 
are subject to the constraint that the volume per 
lipid is conserved. 

With these assumptions, Brannigan and Brown construct 
a free energy functional which is quadratic in the fluctu- 
ations of the fields 2 and A (higher order contributions 
are neglected) and can be used to calculate the thermally 
averaged fluctuations of the total membrane height, h = 
(hi + Ii2)/2, and monolayer thickness, t = (hi — h^jl. 

In our simulations, the situation is different from that 
considered by Brannigan and Brown in two respects. 
First, we apply an external tension. Second, we do not 
have well-separated monolayers, especially at high ten- 
sions. We argue that the second point is not critical: If 
we associate the fields hi_i with the positions of the head 
group layers layers rather than whole monolayers, the as- 
sumptions enumerated above are still reasonable and we 
obtain the same theory. The first point is more subtle, 
since the external tension is not an intrinsic material pa- 
rameter such as, e.g., the bending energy (which drives 
the bending fluctuations). 

To assess the effect of an applied external tension on 
the height fluctuations, we first consider a simplified case, 
where the membrane is characterized by a single surface 
manifold h(x, y) (shifted to (h) = 0) with fixed surface 
area A s and variable projected area A. The surface area 
is related to the projected area via 

A s = J d 2 ? V 1 + (V/0 2 ~^+lJ d 2 r(Vh) 2 . (11) 

Upon applying an external tension T, the free enthalpy 
G of the system is given by 

G = -TA « -TA S + ^ J d 2 r (Vh) 2 , (12) 

where r^4 s is a constant. Hence the external tension 
couples to the fluctuations of the local membrane height 
h(x, y) in the same way as an internal intcrfacial tension 
couples to the fluctuations, e.g., a gas-liquid interface 

The same type of argument can be applied to the model 
of Brannigan and Brown, where the membrane has finite 
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thickness and the lipid volume, rather than the lipid area, 
is conserved: Let z + {x,y) and 2z_{x,y) denote the local 
membrane height and thickness as before, with (z + ) = 0. 
Let furthermore D(x,y) denote the true local membrane 
thickness, evaluated with respect to t he local surface nor- 
mal, i.e., z-(x,y) = D^Jl + (Vz + ) 2 /2. The thickness is 
taken to fluctuate weakly about its mean value D, i.e., 
D(x, y) = D + 5{x, y). The number of lipids on the pro- 
jected area element &A is then given by 



N = 



[ d 2 rz„ = - [ A 2 r Jl + (Vz+) 2 D 

J A V J A 



1 (da + j I 



d 2 r6 



(13) 



with the lipid volume u, where we have expanded up 
to second order in the fluctuating fields h and 6. Upon 
applying an external tension T, the free enthalpy acquires 
an additional term 



G r = -TA = -- 



TvN 



d 2 r(Vz+) 2 + £ 

A L> 



d 2 rS, 
(14) 

where TvN/ D = const. Hence the external tension again 
has the same effect on the height fluctuations than an 
internal tension. The last term in Eq. (|14p results in 
an effective thinning of the membrane and does not con- 
tribute to the fluctuation spectra in the quadratic order 
considered here. 

Supplementing the free energy of Brannigan and 
Brown [36j | with this additional tension term, we obtain 
the Hamiltonian (in Fourier space) 

q 

+ 2(fc A + 7A(? 2 )A+A± q 

+\ E(W*o + Kg 4 - 4fc c C<z 2 AoK *- q 
q 

+ 2(fc A +7A<Z 2 )A q Al q . 

Here 2 ± denotes the bending modes defined above, and 
A ± = (Ai ± A2V2 the corresponding protrusion modes. 
The parameters k c and k A are the bending and compress- 
ibility moduli of the bilayer, £ is related to the sponta- 
neous curvature Co and given by cq — £ dco /d£ with the 
area per lipid S, to is the mean monolayer thickness, and 
the parameters 7^ and k\ characterize the protrusion 
modes. The resulting spectra for height and thickness 
fluctuations are given by 



(15) 



(l^ql 2 ) = 



k R T 



k c q± + Tq 2 2(k x + 7A(7 2 ) 



(16) 



<l*q| 2 } 



k B T 



k B T 



k c q 4 -4k c Cq 2 /t + k A /tl 2{k x +lxq 2 ) 
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FIG. 12: Fluctuation spectra of membranes for different ten- 
sions as indicated. The dashed line is the fit to the thickness 
spectrum. The solid line shows the fit of the elastic theory to 
the height spectrum for fixed tension V and the dotted line 
shows the fit with Ffl t as additional fit parameter. At higher 
tensions these fits deviate from each other. The values of Tfl t 
are listed in Table ITTT1 



We have determined the fluctuation spectra from simu- 
lations of systems containing 3,200 lipids and 24,615 sol- 
vent beads using a method described elsewhere pjl [39[ , 
and used the theory above to fit our data. As in our 
earlier work [16| . the fitting yields very good results for 
tcnsionless membranes. It still works well at the compar- 
atively low tension of T = 1.0 e/ a 2 . In the high-tension 
regime above 2.0e/<7 2 , however, the situation changes. 
When fixing T to the externally applied value, the fit 
consistently underestimates the amplitudes of the long- 
wavelength fluctuations. This can be remedied by leaving 
r as fit parameter. However, the resulting values for the 
fitted tension are smaller by a factor of up to one half 
than the externally applied values (Fig. IT2|) . The fit pa- 
rameters for the fits with fixed and free parameter T are 
given in Tables [XT] and IIIII 

The relation between external (r) and internal stress 
(Fint) in membranes has been discussed by a number of 
authors for the situation where the membrane is kept in 
a frame with fixed projected area A. In these studies, 
the total area A s was allowed to fluctuate, either because 
of fluctuations in the number of molecules (grandcanon- 
ical case) [10, El| or because of fluctuations of the lipid 
area (canonical, compressible case) [12, E3 . The frame 
tension is then found to differ from the internal stress in 
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0.0 



1.0 



k c [e] 


4.9 ± 1.0 


6.2 ±0.5 


c/*o kr 2 ] 


0.12 ±0.02 


0.069 ±0.019 




1.1 ±0.2 


0.92 ±0.09 


kx [e/af] 


1.4 ±0.1 


1.4 ±0.1 


7a [e/a t 2 ] 


0.022 ± 0.039 


0.035 ±0.021 




2.0 


3.0 


kc [e] 


7.0 ±0.3 


7.7 ±0.3 


c/*o kr 2 ] 


-0.038 ± 0.008 


-0.137 ±0.040 


W*o [e/a t 4 ] 


1.8 ±0.2 


4.0 ±0.4 


&A [e/o-f 4 ] 


1.8 ±0.1 


2.0 ±0.1 


7A [e/a t 2 ] 


-0.048 ± 0.042 


-0.062 ± 0.048 



patches. Thus the intrinsic tension should be higher than 
the frame tension, which is opposite from what we ob- 
serve in our simulations. 

However, we believe that the two situations - fixed 
frame and varying surface area versus variable frame and 
(roughly) fixed surface area - are not comparable. Ac- 
cording to the arguments leading to Eqs. (fl~2j) and (fT4| . 
the frame tension is not renormalizcd by fluctuations at 
the level of a Gaussian theory (i.e., a theory based on a 
free energy functional which is quadratic in the fluctua- 
tions). It might be renormalized if one includes higher 
order terms. For example, the last term in Eq. (|14p . 
oc J d 2 r (5, introduces a thickness-mediated interaction 
between the height fluctuation modes via the relation 



S = 2z 



l±(Vz+) 2 -L>«2z_-z_(Vz+) 2 -L>, (18) 



TABLE II: Elastic parameters of the model membrane in the 
fluid phase, obtained from a fit of the elastic theory (Eqs. [TB] 
and EJl with fixed F. 



r [e/af] 


0.0 


1.0 


k c [e] 


4.4 ±0.7 


6.6 ±0.9 


C/*o k" 2 ] 


0.11 ±0.02 


0.074 ± 0.019 


k A /tl [e/at] 


1.0 ±0.2 


0.95 ±0.11 


kx [e/af] 


1.5 ±0.1 


1.4 ±0.1 


7A [e/° 2 t ] 


0.018 ± 0.043 


0.041 ± 0.026 


Ifit [e/a 2 t \ 


0.11 ±0.19 


0.83 ±0.14 




2.0 


3.0 


k c [e] 


10.9 ±0.8 


12.6 ± 1.1 


c/*o kr 2 ] 


-0.011 ±0.012 


-0.11 ±0.08 


k A /tl [e/a?] 


2.0 ±0.3 


4.3 ±0.5 


kx [e/o" t 4 ] 


1.6 ±0.03 


1.7 ±0.1 


7A [e/o- 2 ] 


0.057 ± 0.030 


0.048 ± 0.032 


Tflt [e/o" t 2 ] 


0.90 ±0.31 


1.7 ±0.4 



TABLE III: Elastic parameters of the model membrane in the 
fluid phase from a fit to Eqs. I16l and ll7l where the parameter 
r is allowed to vary. The fitted value is ret. 



the membranes due to the contribution of the membrane 
fluctuations to the surface free energy j40l - {44l |. The cor- 
rection is additive and should always be present, even 
at (external or internal) tension zero. For the canonical, 
compressible case, which is obviously more relevant in 
our context, Farago et al. [4l[ and Imparato 42 1 have 
predicted that the fluctuations reduce the frame tension 
by roughly fc^Tri/A, compared to the intrinsic stress, 
where n is the number of fluctuation degrees of freedom, 
i.e., the number of independently fluctuating membrane 



which might effectively renormalize rj nt . Another pos- 
sibility is of course that the theory of Brannigan and 
Brown j36j . which we have used to analyze the data, is 
no longer applicable at high tensions due to the structural 
rearrangements in the membrane. 



E. Lipid-mediated interactions between inclusions 

Finally in this section, we discuss the effect of ten- 
sion on the membrane-mediated interactions between two 
simple cylindrical inclusions in the bilayer. We focus on 
the effective interactions between these model proteins 
and the influence of an external tension on the potential 
of mean force (PMF). The radial distribution function 
g(r) as a function of the protein-protein distance r was 
obtained from simulation runs using the technique of suc- 
cessive umbrella sampling [4a ] combined with a reweight- 
ing procedure. As starting configurations we used equi- 
librated systems with 750 to 760 lipids and two simple 
transmembrane proteins of diameter 3 a t . A first esti- 
mate of g(r) was obtained during 2 x 10 6 MCS. Then, 
biased runs of 3 x 10 6 MCS were performed to improve 
the statistics of configurationally less frequent protein- 
protein distances. After removing the bias from these 
results and combining the overlapping distributions the 
effective potential w(r) = — ksT In g(r) was extracted. In 
order not to complicate the interpretation of our results, 
the model proteins were not allowed to tilt. This can be 
justified by assuming that real transmembrane proteins 
might be bound to, e.g. cytoskclctal, structures outside 
the membrane, which allow for transverse motion but not 
for tilt. 

The type of model for the inclusions is identical to 
the one introduced in [l6| and a brief overview is given 
in the following: The interaction of this simple model 
protein and the lipid or solvent beads has a repulsive 
contributions, which is described by a radially shifted 
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and truncated Lcnnard- Jones potential 




where r xy = y 1 x 2 + y 2 denotes the distance of the in- 
teraction partners in the xy plane, a is given by a = 
(a t + <7j)/2 for interactions with beads of type i ( i = h, 
t, and s for head, tail, and solvent beads, respectively), 
(To = at, and Vlj has been defined above (Eq.[3]). The di- 
rect protein-protein interactions have the same potential 
(Eq. I19p with a = a t and a$ = 2at- 

In addition, protein cylinders attract tail beads on a 
hydrophobic section of length L. This is described by 
an additional attractive potential that depends on the z 
distance between the tail bead and the protein center. 
The total potential reads 

Vpt( (V Iep (r xy ) + V attl {r X y) X W p (z)) , (20) 

with the attractive Lennard- Jones contribution 

f Vlj(1)-F L j(2) if r xy - a < a 

Kttr = < V L j (l^) - U LJ (2) if a < r xy -a a <2a , 

otherwise 

(21) 

and a weight function W p {z), which is unity on a stretch 
of length 21 = L — 2at and crosses smoothly over to zero 
over a distance of approximately at at both sides. Specif- 
ically, we use 

{1 if \z\ < I 
cos 2 (f|zM) if l< \z\ <; + f . (22) 
otherwise 

The hydrophobicity of the protein is tuned by the pa- 
rameter £ pt . The choice of a sufficiently high interaction 
strength between the hydrophobic core of the membrane 
and the hydrophobic part of the inclusion is crucial to in- 
duce local perturbation of the bilayer. We note that the 
repulsive sections of our model proteins span the whole 
simulation box in z direction. Therefore, the simulations 
were carried out at constant box height L z , volume moves 
were only allowed in the lateral directions x and y, and 
the number of solvent beads of the was allowed to fluc- 
tuate (see Sec. HII- 

We found that the effect of tension on the PMF was 
only significant for rather hydrophobic model proteins, 
i.e., proteins with high interaction parameter e pt . In the 
following, we will present the results obtained with e pt = 
6.0 e. Fig. [T2] shows snapshots of the model proteins in 
membranes at zero tension and at tension T = 2.0 e/a 2 . 

The potentials of mean force between two inclusions 
as a function of their in-plane distance are plotted in 
Fig. [TJ]for two different hydrophobic lengths L. In the 
first case (Fig. [TU top), the inclusion is characterized 
by a negative hydrophobic mismatch in tcnsionless free 




r = 0.0e/o- t 2 r = 2.0e/(T t 2 

FIG. 13: Vertical slices through systems with two inclusions. 
The dark gray rings above the black hydrophobic part of the 
protein indicate the sections where the attractive interaction 
between protein and tail beads (dark gray) decays to zero. 
The light gray section and the light gray beads mark the hy- 
drophilic part of the protein and the lipid heads, respectively. 
The solvent beads above and below the bilayer are marked 
by small dots. The left snapshot shows a stressfree bilayer 
with proteins of hydrophobically matching length L — 6at- 
On the right the same proteins are shown in a bilayer under 
a tension of F = 2.0 e/of. 
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FIG. 14: Potential of mean force of two isolated proteins as 
a function of the protein distance at tension zero (solid lines) 
and r = 2.0 e/of (dotted lines), for proteins with hydropho- 
bic length L — 4<j t (top) and L = 6cr t (bottom). The insets 
show the corresponding radial thickness profiles around a sin- 
gle protein (d p denotes the in-plane distance to the center of 
the inclusion). 



membranes. Under tension, the membrane thins and the 
mismatch is reduced (inset Fig. Q3] top). In the second 
case (Fig. [141 bottom) , the hydrophobic part of the in- 
clusion roughly matches the thickness in the tcnsionless 
case. Under tension, the membrane thins and a posi- 
tive hydrophobic mismatch develops (see inset Fig. 1141 
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bottom). 

Due to lipid packing in the vicinity of the inclusions 
the curves show an oscillatory shape with a wavelength 
of approximately 1 a t . Since at small inclusion- inclusion 
distances direct interactions of the proteins and deple- 
tion induced attraction due to the solvent particles come 
into play, these parts of the curves have been cut off. 
Thus, we can focus on the lipid-mediated medium and 
long ranged interactions. The main effect of tension on 
the lipid mediated interactions between the two model 
proteins can be summarized as follows: In the absence of 
tension, the interactions between hydrophobically mis- 
matched inclusions have an additional attractive contri- 
bution, compared to hydrophobically matched inclusions. 
In the case where the tension reduces the hydrophobic 
mismatch due to membrane thinning, this attraction di- 
minishes. If, on the other hand, the tension leads to a 
stronger hydrophobic mismatch, the attractivity of the 
interaction potential also increases. Therefore, we con- 
clude that the dominant effect of an external tension on 
the lipid-mediated interactions is indirect and related to 
the change of the hydrophobic mismatch due to the mem- 
brane thinning. For decreasing (negative) mismatch, the 
average attraction decreases, and for increasing (positive) 
mismatch, it increases. This is consistent with the behav- 
ior observed in tensionless membranes, where also both 
positive and negative hydrophobic mismatch resulted in 
an attractive contribution to the PMF [l|| . Other recent 
studies [46l - [49| have also highlighted the importance of 
hydrophobic mismatch as a driving factor for protein ag- 
glomeration. From our simulations, an additional effect 
of tension is not evident. 



IV. DISCUSSION AND SUMMARY 

In this paper we have studied the influence of an exter- 
nal tension on the properties of bilayers using a generic 
coarse-grained model. To put these results into perspec- 
tive, we will now briefly discuss the experimental situa- 
tion. 

Experimentally, one of the most widely used techniques 
to determine mechanical stretch properties of bilayers 
is the micropipctte approach, where giant bilayer vesi- 
cles are pressurized by micropipette suction [50j. This 
method produces a uniform membrane tension, is very 
accurate, and can be used to verify elastic reversibil- 
ity [5l[. Micropipette aspiration experiments, e.g., car- 
ried out by Ncedham and Nun [52j, found for different 
lipids and lipid/cholcsterol mixtures that membrane lysis 
is usually reached at an relative areal expansion of less 
than 5 %, and the rupture strength at low cholesterol 
concentration was typically around 2 — lOmN/m. 

This does not compare well with our findings and those 
of other atomistic or mesoscopic simulations (see the in- 
troduction to Sec. 1111)1 . where fluid membranes could sus- 
tain tensions of 30mN/m or more and remained stable 
up to relative extensions of 40 % or more. It should how- 



ever be noted that tension-induced lysis is a stochastic 
process, and the duration of the exposure to the stress 
plays a decisive role for the stability. Experimentally, 
abrupt failure of the bilayer is observed, when the ten- 
sion reaches a critical value, whereas below this tension 
long-term persistence of the stressed membrane can be 
presumed [53j. Evans et al. showed that the rupture 
strength of membranes is a p roperty which crucially de- 
pends on the loading rate [54j . In their study on five types 
of fluid giant phosphatidylcholine lipid vesicles, they var- 
ied the loading rate from 0.01 — lOOmN/m/s. At high 
loading rates the systems were found to be stable up to 
values of 20 - 30mN/m. 

To get a rough estimate of the meaning of the 
timescales in our simulations, compared to experimen- 
tal systems, we can map the diffusion constants D xy 
in Sec. IIII CI to the corresponding values measured for 
DPPC bilayers in recent experiments. At 45° C Scom- 
parin et al. [55| report a diffusion coefficient of ap- 
proximately 10/im 2 /s. Taking our result of D xy = 
0.8 x 10~ 5 of /MCS for the tensionless membrane and set- 
ting our intrinsic length scale 1 at to 6 A as described ear- 
lier, we find that 10 6 MCS in our simulation correspond 
to approximately 0.3 (is in real time. The typical lengths 
of our simulations lie between 2-4 x 10 6 MCS, which cor- 
responds to roughly 1 fxs in real systems. Thus our simu- 
lations correspond to systems exposed to very high load- 
ing rates and short timescales. Taking this into consid- 
eration, their stability does not contradict experimental 
findings. Our simulations provide a way to study mem- 
branes under extreme conditions, and to analyze their 
structural properties, which cannot be accessed easily by 
experiments. 

Due to these difficulties, experimental results with 
which we could compare our simulations are scarce. One 
positive example is the behavior of the the ripple phase 
under stress. To our knowledge, we have performed the 
first simulation study which tries to shed light on the 
structural rearrangements and elastic properties of a bi- 
layer in the ripple phase under lateral stress. We have 
shown that lateral tension leads to suppression of the rip- 
ple structure in the Pp* phase, and we see a transition of 
the areal extensibility from soft, fluid-like to gel-like be- 
havior. Qualitatively, this behavior agrees well with the 
findings of Evans et al. [53|, |j56[ , who also reported an ini- 
tial soft-elastic response at low tensions and stiff clastic 
properties after elimination of the ripple. 

In the gel phase, the tension does not change the state 
of the bilayer significantly in the range of tensions consid- 
ered in this work. The basic structure of the two mono- 
layers stays intact. The situation is very different in the 
biologically most relevant fluid phase. At temperatures 
where the tensionless membranes are fluid, they respond 
to high tensions by a structural change from a state where 
both monolayers are well-separated to a state where they 
are partly interdigitatcd. These changes are associated 
with substantial variations of the compressibility (up to 
a factor of 3), and the lipid diffusion constant (up to 20 
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%). 

We have also studied the influence of membrane ten- 
sion on the effective interaction between two model pro- 
teins. Under tension, the membrane becomes thin- 
ner, which affects the hydrophobic mismatch interaction. 
This is found to be the dominant effect. The interaction 
between negatively mismatched proteins decreases, and 
that between positively mismatched proteins increases. 
Thus applying tension can be used to tune the strength 
of membrane mediated protein-protein interactions. 
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